6 results
Analysis of scale-dependent kinetic and potential energy in sheared, stably stratified turbulence
- Xiaolong Zhang, Rohit Dhariwal, Gavin Portwood, Stephen M. de Bruyn Kops, Andrew D. Bragg
-
- Journal:
- Journal of Fluid Mechanics / Volume 946 / 10 September 2022
- Published online by Cambridge University Press:
- 29 July 2022, A6
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
-
Budgets of turbulent kinetic energy (TKE) and turbulent potential energy (TPE) at different scales $\ell$ in sheared, stably stratified turbulence are analysed using a filtering approach. Competing effects in the flow are considered, along with the physical mechanisms governing the energy fluxes between scales, and the budgets are used to analyse data from direct numerical simulation at buoyancy Reynolds number $Re_b=O(100)$. The mean TKE exceeds the TPE by an order of magnitude at the large scales, with the difference reducing as $\ell$ is decreased. At larger scales, buoyancy is never observed to be positive, with buoyancy always converting TKE to TPE. As $\ell$ is decreased, the probability of locally convecting regions increases, though it remains small at scales down to the Ozmidov scale. The TKE and TPE fluxes between scales are both downscale on average, and their instantaneous values are correlated positively, but not strongly so, and this occurs due to the different physical mechanisms that govern these fluxes. Moreover, the contributions to these fluxes arising from the sub-grid fields are shown to be significant, in addition to the filtered scale contributions associated with the processes of strain self-amplification, vortex stretching and density gradient amplification. Probability density functions (PDFs) of the $Q,R$ invariants of the filtered velocity gradient are considered and show that as $\ell$ increases, the sheared-drop shape of the PDF becomes less pronounced and the PDF becomes more symmetric about $R=0$.
Hydrodynamic interactions and extreme particle clustering in turbulence
- Andrew D. Bragg, Adam L. Hammond, Rohit Dhariwal, Hui Meng
-
- Journal:
- Journal of Fluid Mechanics / Volume 933 / 25 February 2022
- Published online by Cambridge University Press:
- 23 December 2021, A31
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
-
Expanding recent observations by Hammond & Meng (J. Fluid Mech., vol. 921, 2021, A16), we present a range of detailed experimental data of the radial distribution function (r.d.f.) of inertial particles in isotropic turbulence for different Stokes number, $St$, showing that the r.d.f. grows explosively with decreasing separation r, exhibiting $r^{-6}$ scaling as the collision radius is approached, regardless of $St$ or particle radius $a$. To understand such explosive clustering, we correct a number of errors in the theory by Yavuz et al. (Phys. Rev. Lett., vol. 120, 2018, 244504) based on hydrodynamic interactions between pairs of small, weakly inertial particles. A comparison between the corrected theory and the experiment shows that the theory by Yavuz et al. underpredicts the r.d.f. by orders of magnitude. To explain this discrepancy, we explore several alternative mechanisms for this discrepancy that were not included in the theory and show that none of them are likely the explanation. This suggests new, yet-to-be-identified physical mechanisms are at play, requiring further investigation and new theories.
Clustering of rapidly settling, low-inertia particle pairs in isotropic turbulence. Part 2. Comparison of theory and DNS
- Sarma L. Rani, Rohit Dhariwal, Donald L. Koch
-
- Journal:
- Journal of Fluid Mechanics / Volume 871 / 25 July 2019
- Published online by Cambridge University Press:
- 22 May 2019, pp. 477-488
-
- Article
- Export citation
-
Part 1 (Rani et al. J. Fluid Mech., vol. 871, 2019, pp. 450–476) of this study presented a stochastic theory for the clustering of monodisperse, rapidly settling, low-Stokes-number particle pairs in homogeneous isotropic turbulence. The theory involved the development of closure approximations for the drift and diffusion fluxes in the probability density function (p.d.f.) equation for the pair relative positions $\boldsymbol{r}$. In this part 2 paper, the theory is quantitatively analysed by comparing its predictions of particle clustering with data from direct numerical simulations (DNS) of isotropic turbulence containing particles settling under gravity. The simulations were performed at a Taylor micro-scale Reynolds number $Re_{\unicode[STIX]{x1D706}}=77.76$ for three Froude numbers $Fr=\infty ,0.052,0.006$, where $Fr$ is the ratio of the Kolmogorov scale of acceleration and the magnitude of gravitational acceleration. Thus, $Fr=\infty$ corresponds to zero gravity, and $Fr=0.006$ to the highest magnitude of gravity among the three DNS cases. For each $Fr$, particles of Stokes numbers in the range $0.01\leqslant St_{\unicode[STIX]{x1D702}}\leqslant 0.2$ were tracked in the DNS, and particle clustering quantified both as a function of separation and the spherical polar angle. We compared the DNS and theory values for the exponent $\unicode[STIX]{x1D6FD}$ characterizing the power-law dependence of clustering on separation. The $\unicode[STIX]{x1D6FD}$ from the $Fr=0.006$ DNS case are in reasonable agreement with the theoretical predictions obtained using the second drift closure (referred to as DF2). To quantify the anisotropy in clustering, we calculated the leading–order coefficient in the spherical harmonics expansion of the p.d.f. of pair relative positions. The coefficients predicted by the theory (DF2) again show reasonable agreement with those calculated from the DNS clustering data for $Fr=0.006$. However, we note that in spite of the high magnitude of gravity, the clustering is only marginally anisotropic both in DNS and theory. The theory predicts that the spherical harmonic coefficient scales with $\unicode[STIX]{x1D6FD}(=\unicode[STIX]{x1D6FD}_{2}St_{\unicode[STIX]{x1D702}}^{2})$, where $\unicode[STIX]{x1D6FD}_{2}$ is the ratio of the drift and diffusion flux coefficients. Since the drift flux, and thereby $\unicode[STIX]{x1D6FD}_{2}$, is seen to decrease with gravity for $St_{\unicode[STIX]{x1D702}}<1$, the anisotropy is also correspondingly diminished.
Small-scale dynamics of settling, bidisperse particles in turbulence
- Rohit Dhariwal, Andrew D. Bragg
-
- Journal:
- Journal of Fluid Mechanics / Volume 839 / 25 March 2018
- Published online by Cambridge University Press:
- 02 February 2018, pp. 594-620
-
- Article
- Export citation
-
Mixing and collisions of inertial particles at the small scales of turbulence can be investigated by considering how pairs of particles move relative to each other. In real problems the two particles will have different sizes, i.e. they are bidisperse, and the effect of gravity on their motion is often important. However, how turbulence and gravity compete to control the motion of bidisperse inertial particles is poorly understood. Motivated by this, we use direct numerical simulations (DNS) to investigate the dynamics of settling, bidisperse particles in isotropic turbulence. In agreement with previous studies, we find that without gravity (i.e. $Fr=\infty$, where $Fr$ is the Froude number), bidispersity leads to an enhancement of the relative velocities, and a suppression of their spatial clustering. For $Fr<1$, the relative velocities in the direction of gravity are enhanced by the differential settling velocities of the bidisperse particles, as expected. However, we also find that gravity can strongly enhance the relative velocities in the ‘horizontal’ directions (i.e. in the plane normal to gravity). This non-trivial behaviour occurs because fast settling particles experience rapid fluctuations in the fluid velocity field along their trajectory, leading to enhanced particle accelerations and relative velocities. Indeed, the results show that even when $Fr\ll 1$, turbulence can still play an important role, not only on the horizontal motion, but also on the vertical motion of the particles. This is related to the fact that $Fr$ only characterizes the importance of gravity compared with some typical acceleration of the fluid, yet accelerations in turbulence are highly intermittent. As a consequence, there is a significant probability for particles to be in regions of the flow where the Froude number based on the local, instantaneous fluid acceleration is ${>}1$, even though the typically defined Froude number is $\ll 1$. This could imply, for example, that extreme events in the mixing of settling, bidisperse particles are only weakly affected by gravity even when $Fr\ll 1$. We also find that gravity drastically reduces the clustering of bidisperse particles. These results are strikingly different to the monodisperse case, for which recent results have shown that when $Fr<1$, gravity strongly suppresses the relative velocities in all directions, and can enhance clustering. Finally, we consider the implications of these results for the collision rates of settling, bidisperse particles in turbulence. We find that for $Fr=0.052$, the collision kernel is almost perfectly predicted by the collision kernel for bidisperse particles settling in quiescent flow, such that the effect of turbulence may be ignored. However, for $Fr=0.3$, turbulence plays an important role, and the collisions are only dominated by gravitational settling when the difference in the particle Stokes numbers is ${\geqslant}O(1)$.
Stochastic theory and direct numerical simulations of the relative motion of high-inertia particle pairs in isotropic turbulence
- Rohit Dhariwal, Sarma L. Rani, Donald L. Koch
-
- Journal:
- Journal of Fluid Mechanics / Volume 813 / 25 February 2017
- Published online by Cambridge University Press:
- 17 January 2017, pp. 205-249
-
- Article
- Export citation
-
The relative velocities and positions of monodisperse high-inertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al., J. Fluid Mech., vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian two-time correlation of fluid relative velocities seen by pairs that are nearly stationary. The two-time correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by large-scale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integral-scale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS. LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokker–Planck form of the pair PDF equation. Langevin simulations of particle-pair dispersion were performed using three closure forms of the diffusivity – i.e. the one containing the time integral of the Eulerian two-time correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the two-time correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relative-motion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particle-laden forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$. Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$ is the Taylor micro-scale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particle-pair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(U|r)$ and $\unicode[STIX]{x1D6FA}(U_{r}|r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$. The first closure form, involving the Eulerian two-time correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs.
A stochastic model for the relative motion of high Stokes number particles in isotropic turbulence
- Sarma L. Rani, Rohit Dhariwal, Donald L. Koch
-
- Journal:
- Journal of Fluid Mechanics / Volume 756 / 10 October 2014
- Published online by Cambridge University Press:
- 05 September 2014, pp. 870-902
-
- Article
- Export citation
-
The probability density function (PDF) kinetic equation describing the relative motion of inertial particle pairs in a turbulent flow requires closure of the phase-space diffusion current. A novel analytical closure for the diffusion current is presented that is applicable to high-inertia particle pairs with Stokes numbers $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}{\mathit{St}}_r \gg 1$. Here ${\mathit{St}}_r$ is a Stokes number based on the time scale $\tau _r$ of eddies whose size scales with pair separation $r$. In the asymptotic limit of ${\mathit{St}}_r \gg 1$, the pair PDF kinetic equation reduces to an equation of Fokker–Planck form. The diffusion tensor characterizing the diffusion current in the Fokker–Planck equation is equal to $1/\tau _v^2$ multiplied by the time integral of the Lagrangian correlation of fluid relative velocities along particle-pair trajectories. Here, $\tau _v$ is the particle viscous relaxation time. Closure of the diffusion tensor is achieved by converting the Lagrangian correlations of fluid relative velocities ‘seen’ by pairs into Eulerian fluid-velocity correlations at pair separations that remain essentially constant during time scales of $O(\tau _r)$; the pair centre of mass, however, is not stationary and responds to eddies with time scales comparable to or smaller than $\tau _v$. For isotropic turbulence, Eulerian fluid-velocity correlations may be expressed as Fourier transforms of the velocity spectrum tensor, enabling us to derive a closed-form expression for the diffusion tensor. A salient feature of this closure is that it has a single, unique form for pair separations spanning the entire spectrum of turbulence scales, unlike previous closures that involve velocity structure functions with different forms for the integral, inertial subrange, and Kolmogorov-scale separations. Using this closure, Langevin equations, which are statistically equivalent to the Fokker–Planck equation, were solved to evolve particle-pair relative velocities and separations in stationary isotropic turbulence. The Langevin equation approach enables the simulation of the full PDF of pair relative motion, instead of only the first few moments of the PDF as is the case in a moments-based approach. Accordingly, PDFs $\varOmega (U|r)$ and $\varOmega (U_r|r)$ are computed and presented for various separations $r$, where the former is the PDF of relative velocity $U$ and the latter is the PDF of the radial component of relative velocity $U_r$, both conditioned upon the separation $r$. Consistent with the direct numerical simulation (DNS) study of Sundaram & Collins (J. Fluid Mech., vol. 335, 1997, pp. 75–109), the Langevin simulations capture the transition of $\varOmega (U|r)$ from being Gaussian at integral-scale separations to an exponential PDF at Kolmogorov-scale separations. The radial distribution functions (RDFs) computed from these simulations also show reasonable quantitative agreement with those from the DNS study of Février, Simonin & Legendre (Proceedings of the Fourth International Conference on Multiphase Flow, New Orleans, 2001).